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Abstract 

In  this  work,  cold  start-up  of  a  methane  fuel  processor  is  explored.  The  experimental  fuel  processor  is  intended  to  provide  hydrogen  for  a 
proton  exchange  membrane  (PEM)  fuel  cell  for  the  power  generation  (3  kWe).  A  dynamic  model  describing  a  series  of  reactors,  the  reformer,  three 
water-gas  shift  reactors,  and  preferential  reactor  is  constructed.  Two  important  factors  for  rapid  start-up  are  identified:  speed  of  temperature  front 
propagation  and  acceptable  CO  concentration.  Steady- state  analyses  reveal  that  the  fuel  feed  flow  rate  with  fixed  steam- to-carbon  and  air-to-carbon 
ratios  is  an  ideal  manipulated  variable.  Considering  both  large  initial  heat  flux  and  gradual  transition  back  to  nominal  operation,  the  shape  of 
feed  manipulation  is  determined.  With  the  feed  scenario  available,  the  fuel  processor  start-up  can  be  formulated  as  a  constrained  optimization 
problem  and  can  be  solved  numerically.  From  optimization  result,  a  heuristic  is  generated  for  rapid  start-up  of  a  fuel  processor.  This  leads  to  a  25% 
improvement  in  the  start-up  time.  Finally,  issues  of  design  modification  are  explored  for  further  reduction  in  the  start-up  time. 

©  2006  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Fuel  cell  systems  offer  high  potential  for  efficiency  and 
reduced  emissions  in  power  generation  [1].  The  proton  exchange 
membrane  fuel  cell  (PEMFC)  is  one  of  the  most  popular  fuel  cell 
systems  in  which  fuels  such  as  methanol  or  methane  is  converted 
to  hydrogen  rich  syn-gas  in  a  reformer  and  which  is  subsequently 
used  in  the  fuel  cell  stack.  In  addition  to  the  reformer,  a  series 
of  CO  reducing  steps,  water  gas  shift  reactions  and  preferential 
oxidation  reactions  were  taken  to  keep  CO  concentration  below 
100  ppm  before  the  syn-gas  entering  the  cell  stack.  This  combi¬ 
nation  constitutes  the  entire  fuel  processor  [2] .  A  dynamic  model 
is  essential  for  the  fuel  processor  operation  for  the  following  rea¬ 
sons:  (1)  discriminating  control  system  design  for  improved  load 
rejection,  and  (2)  evaluating  start-up  strategies  for  fast  start-up. 

Extensive  literature  has  examined  various  aspects  of  fuel 
processor  systems  for  hydrogen  rich  syn-gas  production  which 
include  overviews  of  the  fuel  processing  technology  [3-8],  in 
which  the  reforming  technology  of  hydrocarbon  fuels  is  still  the 
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major  focus.  Steady-state  simulations  are  often  performed  for 
sensitivity  analyses  in  the  design  and  operation  phases  of  the 
fuel  processor  [9-12].  Studies  on  dynamic  behavior  of  the  fuel 
processor  have  received  some  attention  lately,  and  typically  the 
relationships  between  feed  conditions  and  dynamic  responses 
were  explored  in  [13,14].  The  start-up  dynamics  was  explored 
in  [15]  in  order  to  devise  a  more  efficient  start-up  strategy. 

The  objective  of  this  work  is  to  construct  a  dynamic  model 
for  a  methane  fuel  processor  and  control  structure  can  evaluate 
the  performance  at  start-up  stage.  The  remainder  of  this  paper  is 
organized  as  follows.  Section  2  describes  process  and  modeling 
aspects  of  the  fuel  processor.  Sensitivities  of  operating  parame¬ 
ters  are  explored  in  Section  3,  followed  by  start-up  strategy  and 
evaluation  in  Section  4.  Conclusion  is  drawn  in  Section  5. 

2.  Process  studied 

An  experimental  methane  fuel  processor  is  constructed  to 
generate  hydrogen  for  3  kWe  fuel  cell  application.  Fig.  1  shows 
the  schematic  of  the  fuel  processor  system  with  corresponding 
dimensions  (marked  in  number  with  the  unit  of  cm).  The  reac¬ 
tants,  methane,  air,  and  water,  pass  through  two  heat  exchangers 
right  after  entering  the  system  as  shown  in  the  upper  left  corner 
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Nomenclature 

Cp  heat  capacity  of  the  gas  (kJ  mol  K- 1 ) 

Cp7 s  heat  capacity  of  the  carrier  (kJ  mol  K- 1 ) 

Cpw  heat  capacity  of  the  metal  reforming  reactor  wall 
(kJmol  K-1) 

D\  inner-diameter  of  the  reactor  (cm) 

D0  outer-diameter  of  the  reforming  reactor  (cm) 

F  total  molar  flow  rate  (mol  min- 1 ) 
h  pulse  height 

A77R  heat  of  reaction  of  reaction  (kJ  mol-1) 
kCond  thermal  conductivity  of  the  metal  reactor  wall 
(kJ  min-1  cm-1  K-1) 

mw  weight  of  the  reforming  reactor  (g) 

Mp  molar  holdup  of  the  burner  (mol) 

P  pressure  (atm) 

<2f  heat  input  for  preheating  (kJ  min- 1 ) 

r  rate  of  reaction  (mol  g- 1  min- 1 ) 

S  ramp  down  slope 

td  pulse  duration  time 

T  reaction  temperature  (K) 

7a  surrounding  temperature  (K) 

7f  temperature  of  the  feed  (K) 

7jn  inlet  temperature  of  the  reformer  (K) 

Thi  inlet  temperature  of  the  HTS1  (K) 

7h2  inlet  temperature  of  the  HTS2  (K) 

7l  inlet  temperature  of  the  LTS  (K) 

7p  inlet  temperature  of  the  PROX  (K) 

7w  reactor  wall  temperature  (K) 

U  heat  transfer  coefficient  (kJ  min- 1  cm-2  K- 1 ) 

Vr  volume  of  the  gas  (cm3) 

Wcat  catalyst  weight  (g) 

Ws  carrier  weight  (g) 

y  mole  fraction 

Greek  letters 

v  stoichiometric  coefficient  of  the  reaction 

p  density  of  the  carrier  (g  cm-3) 

pav  average  density  of  the  gas  in  the  reforming  reactor 
(g  cm-3) 


of  Fig.  1 .  Next,  the  mixed  reactant  goes  through  the  coil  and 
then  the  burner  before  entering  the  ATR.  The  ATR  effluent  is 
passed  to  the  gas  cleaning  units  (HTS 1 ,  HTS2,  PROX1-PROX3) 
as  shown  in  Fig.  1.  In  the  modeling  phase,  the  experimental 
fuel  processor  is  simplified  to  a  reformer,  a  burner,  three  water 
gas  shift  reactors  (HTS1,  HTS2,  LTS),  and  a  preferential  oxida¬ 
tion  reactor  (PROX)  as  shown  in  Fig.  2.  Fresh  feed  containing 
methane,  air,  and  water  is  fed  into  the  reformer  to  carry  out 
autothermal  reforming  (ATR).  The  reformer  was  integrated  with 
a  burner  which  supplies  the  heat  needed  for  the  reactions.  The 
reformate  is  passed  through  a  feed-effluent  heat  exchanger  to 
the  first  high-temperature  water-gas  shift  reactor  (HTS1)  fol¬ 
lowed  by  the  second  reactor  (HTS 2)  and  the  to  a  low- temperature 
water-gas  shift  reactor  (LTS).  In  these  three  reactors,  CO  was 


removed  to  meet  the  specification.  Because  the  CO  concentra¬ 
tion  out  of  the  LTS  was  still  too  high,  the  preferential  oxidation 
reaction  (PROX)  is  performed.  In  the  PROX  reactors  O2  injec¬ 
tion  devices  are  installed  at  the  inlet  of  PROX,  and  CO  is  further 
oxidized  to  CO2,  while,  at  the  same  time,  some  H2  is  converted 
to  H2O. 

2.7.  Reaction  kinetics 

Table  1  shows  the  reactions  occurred  in  the  reformer,  HTS1 
and  HTS2,  LTS,  and  PROX.  As  suggested  by  Lin  et  al.  [16], 
a  steady- state  model  describing  the  series  of  reactors  has  been 
constructed.  With  available  experimental  inlet  and  outlet  tem¬ 
peratures  and  compositions,  the  kinetic  parameters  have  been 
adjusted  to  describe  experimental  data.  Let  us  use  the  reformer 
as  an  example.  The  pre-exponential  factors  are  adjusted  to  match 
the  experimental  results  while  keeping  the  activation  energies 
and  adsorption  constants  unchanged.  Hence,  the  data  fitting 
algorithm  of  the  reformer  has  the  following  steps: 

(1)  Starting  with  the  SR  pathway,  the  reactor  wall  inlet  temper¬ 
ature,  7w,in>  is  changed  to  match  the  outlet  temperature  of 
SR  and  the  pre-exponential  factors  of  r\  and  r2  are  adjusted 
to  match  the  outlet  component  flow  rates. 

(2)  While  keeping  the  kinetics  parameter  of  r\  and  r2  fixed, 
7w,in  and  pre-exponential  factor  of  r 3  are  varied  to  match 
the  outlet  temperature  and  component  flows  for  the  ATR 
pathway. 

(3)  Constrain  7w,in  by  establishing  a  linear  equation  to  relate 
7w,in(0)  with  Tm. 

(4)  Modify  kinetics  parameters  in  r\  and  r3  such  that  the  errors 
are  of  the  same  order  of  magnitude  for  the  SR  and  ATR 
pathways. 

The  procedure  is  repeated  for  HTS1,  HTS2,  LTS  and  PROX 
reactors.  The  rate  expressions  and  corresponding  parameter  val¬ 
ues  are  shown  in  Table  1  for  the  entire  fuel  processor. 

2.2.  Dynamic  modeling 

Homogenous  reactor  models  were  set  up  to  describe  the 
dynamic  behavior  of  the  experimental  fuel  processor  [17].  The 
assumptions  made  include: 

1 .  Constant  pressure  ( 1  bar). 

2.  Ideal  gas  behavior. 

3.  Same  temperature  for  the  vapor  and  solid  phases. 

4.  Negligible  thermal  capacity  for  the  gas  phase  in  the  reactor 
as  compared  to  that  of  the  solid  catalyst  and  carrier. 

5.  Negligible  heat  loss  for  the  HTS1,  HTS2,  LTS  and  PROX 
reactors. 

Partial  differential  equations  describing  the  energy  and  mass 
balances  are  lumped  into  N  sections  using  the  honeycomb  car¬ 
rier  weight  (Ws)  as  the  independent  variable  (note  that  this  also 
applies  to  the  catalyst  weight  Wcat,  become  it  is  assumed  that  the 
catalyst  is  distributed  uniformly  throughout  the  carrier).  Thus, 
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Fig.  1.  Experimental  fuel  processor  and  corresponding  dimensions. 


Fig.  2.  Process  configuration  and  control  structure  of  the  fuel  processor. 
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Table  1 


Reaction  rate  expressions  and  parameter  values  for  the  fuel  processor  system 


Reactor 

Reaction 

Kinetics 

AHr  (kJ  mol  1 ) 

ATR 

CH4  +  H20  CO  +  3H2 

k\  Pcn4  Pn2o  ~  k\  Pqo P^ 

n  ~ 

206 

C0  +  H20  4>C02  +  H2 

h  pCO  Pw  0~k'2  Pc02  Pro 

n  -  Ph2<“>>2 

-41.2 

CH4  +  202  =>•  C02  +  2H20 

_  ha  Pcn4  Po2  hbPcn4Po2 

r3  ~~  («2)2  +  «2 

-810 

HTS1,  HTS2 

C0  +  H20  4»C02  +  H2 

m ts  =  Pqo  Ph2o  —  kfH  P\\2  Pco2 

-41.2 

LTS 

co+h2o^co2+h2 

0_ts  =  ^lRco^h2o  -  k'HPn2Pco2 

-41.2 

PROX 

co  +  \o2  =>  co2 

^PROXl  =  kcoPcO 

-283 

h2  +  2o2  =:>  h2o 

rPROX2  =  1.5rco 

-243 

(xi  =  ^1  +  KqoPco  +  Kh2  pn2  +  Kch4  Pcn4  +  Kh2o  >a  2  =  (1  +  ^ch4  ^ch4  +  Kofpo2)  with  Kqh4  =  6.65  x  10  4  e4607/r,  Kqo  =  8.23  x  10  5  e8504/r, 
KHl  =  6.12  x  10“9  e"71/r,  KH2q  =  1.77  x  105  e-10669/r,  K^4  =  1.26  x  10_1  e3284/r,  and  =  7.87  x  10-7  e11162/r. 


we  have  N  ordinary  differential  equations  instead  of  one  par¬ 
tial  differential  equation  for  each  component  material  balance 
[18].  Consider  the  nth  section  in  the  axial  direction.  The  energy 
balance  of  the  reformer  can  be  expressed  as: 


the  thermal  conductivity  of  the  reactor  wall,  and  D0  is  the  outer- 
diameter  of  the  reforming  reactor.  7a  represents  the  ambient 
temperature.  The  inlet  temperature  of  the  reformer  is  heated  by 
a  burner  and  the  equation  becomes 


,  —  Fn—\Cp,n—lTn—\  FnC p,nTn 


d  t 


-^»EAflR.Ar 


4 U(Tn  -  Tw,„)Ws,„ 


1 


PsDi 


(1) 


where  CPi s  is  the  heat  capacity  of  the  carrier,  Ws,n  the  weight 
of  the  carrier,  Tn  and  7w,/?  the  reaction  temperature  and  reactor 
wall  temperature  in  the  nth  lump,  respectively,  Fn  the  total  molar 
flow  rate  at  the  nth  lump,  Cp?n  the  heat  capacity  of  the  gas  in  the 
nth  lump,  Wcat./i  the  weight  of  the  catalyst  in  the  nth  lump,  rnj 
the  reaction  rate  of  the  jth  reaction  at  the  nth  section,  A77rj  the 
heat  of  reaction  for  the  y'th  reaction,  U  the  overall  heat  transfer 
coefficient,  D\  the  inner-diameter  of  the  reactor,  and  ps  is  the 
density  of  the  carrier.  The  component  material  balance  for  the 
ith  composition  becomes 


dyn  i  v — 

Pav^R,n~ =  Fn—\yn—\,i  Fnyn,i  T  f^cat ,n  /  (2) 


where  pav  is  the  averaged  density  of  the  gas  in  the  reforming 
reactor,  Vr^  the  volume  of  the  gas  in  the  nth  lump  of  the  reformer, 
i>ij  the  stoichiometric  coefficient  of  the  ith  component  under  jth 
reaction,  and  ynj  is  the  mole  fraction  of  the  ith  component  in  the 
nth  lump.  Similarly,  the  reactor  metal  wall  temperature  can  be 
expressed  as 


tnw.nCpWji 


d7w, 


n 


d  t 


=  fccond^L(7w,w-l  ?W  ,n) 


+ 


+ 


4U(Tn  -  Tw  n)Ws, 


n 


pD\ 

4£/(7a  -  Tw,n)Ws,n 


pD 


o 


where  mw,n  is  the  weight  of  the  metal  reactor  wall  in  the  nth 
lump,  Cpw,«  the  heat  capacity  of  the  metal  reactor  wall,  kcond 


Cp.avgMp^  =  Q¥  +  fY/  y,Cp.,(Tt  -  Tin)  (4) 

where  CpAYg  is  the  averaged  heat  capacity  of  the  feed,  Mp  the 
molar  holdup  of  the  burner,  gp  the  heat  needed  for  preheating, 
and  7f  is  the  temperature  of  the  fresh  feed.  Relationship  between 
the  reactor  inlet  temperature  T4n  (To  is  the  lumped  notation)  and 
the  reactor  wall  temperature  at  the  inlet  (7w,o  is  the  lumped 
notation)  is  established  from  a  regression  model  of  [16].  The 
energy  balance  equations  describing  the  burner  provide  the  inlet 
conditions  for  the  reformer  gas  and  metal  wall  temperatures. 

Similarly,  the  equations  describing  water  gas  shift  reactor  and 
PROX  can  be  derived.  The  energy  balance  equation  becomes 


Cp,sWs,n 


d  T 
d  t 


—  Fn  —  1  Cp^n  —  1  Tn— 1  FnCp,nTn 

—  Wcat,  n  A  Hr  j  Vjj  rn ,  j 


(5) 


and  the  component  material  balance  equation  is 

dyn ,  i  \ — 'v 

Pav^R,n~ —  Fn  —  iyn  —  l,i  Fnyn,i  T  ^cat,n/  j^ijrn,j  (6) 


1 


The  modeling  equations  of  the  HTS1,  HTS2,  LTS  and  PROX 
were  assumed  to  be  adiabatic  and  they  were  simpler  than  the 
modeling  equation  of  the  reformer. 

The  rate  expressions  of  the  reactions  (rj)  that  take  place  in  the 
fuel  processor  are  obtained  form  the  regression  of  the  experimen¬ 
tal  data  [14].  For  the  sake  of  completeness,  the  rate  expressions 
are  given  in  Table  1  and  corresponding  parameter  values  are 
shown  in  Table  2.  Table  3  summarizes  steady- state  operating 
condition  with  H2O/CH4  feed  ratio  of  1.45,  O2/CH4  feed  ratio 
of  0.45  while  the  reformer  inlet  temperature  was  set  to  717  °C 
[14]. 

Before  getting  into  the  simulation  results,  it  is  also  impor¬ 
tant  to  note  that  a  heat-integration  scheme  can  be  applied  to 
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Table  2 


Kinetic  data  from  regression  of  an  experimental  methane  fuel  processor 


Reactor 

Reaction 

Pre-exponential 
factor,  <2o 

Activation  energy 
(kJ  mol- 1 ) 

ATR 

r(l),  forward  k\ 

6.32  x  1016 

240.1 

r(l),  reverse  k\ 

1.759  x  103 

17.0 

r( 2),  forward 

2.77  x  106 

67.1 

r(3),  forward 

1.56  x  10s 

86.0 

r(3),  forward  k^ 

1.31  x  10s 

86.0 

HTS 

rHTS  >  forward  k\\ 

9.886  x  105 

47.4 

rHTS,  reverse  k'H 

1.32  x  10“2 

38.1 

LTS 

/"lts  ,  forward  &l 

1.285  x  106 

47.4 

0_ts  » reverse  k'L 

1.32  x  10“2 

38.1 

PROX 

rpROx,  forward  &prox 

1.34  x  104 

8.3 

O 

o 


700 

600 

500 

400 

300 

200 

100 

0 


ATR 

HTS1 

HTS2 

LTS 

PROX 

c 

“ 

\/ 

\  v 

\  \ 
H1° 

V 

■*  i 

*  • 

*  » 

»  % 

■  i  t 

l  ' 

A 

Is; 

t=teo 

■ 

• 

A 

« 

A 

t=45 

-—-cm 

% 

% 

% 

\  t 

t=25  » 

IA30 

\  ...  % 

»  • 

•  V  * 

1  \  * 
t=0  njiin  \  \ 

•  i 

t  i 

%  i 

*  « 

\  % 

% 

%  : 

_ 

_ _ ! _ ! _ 1 _ ! _ 

0 


(A) 


0.2  0.4 

Length  (m) 


0.6 


0.8 


save  the  fuel  consumption.  This  is  especially  true  at  the  start-up 
stage.  The  reason  is  that  before  the  CO  concentration  reaches 
45  ppm  level,  hydrogen  is  recycled  back  to  the  burner  to  be  used 
as  the  fuel  while  the  remaining  portion  of  hydrogen  is  used  to 
preheat  the  feed.  Fig.  1  shows  that  the  inlet  temperature  to  the 
reformer  is  controlled  by  manipulating  the  fuel  flow  rate  while 
a  valve  position  control  is  installed  to  ensure  minimal  utiliza¬ 
tion  of  external  fuel  by  adjusting  the  hydrogen  flow  rate  into  the 
burner  while  the  remaining  hydrogen  is  passed  through  feed- 
effluent  heat  exchanger  (FEHE)  to  provide  additional  preheating 
[19]. 

2.3.  Simulation  results 

The  ordinary  differential  equations  (ODE;  Eqs.  ( 1)— (6)) 
describing  the  fuel  processor  are  solved  as  the  initial  value 
problem  using  Euler  method.  The  ODEs  are  programmed  in 
FORTRAN  and  a  typical  simulation  run  takes  about  2  min  using 
IBM  PC  with  Pentium  4  processor. 


Table  3 

Steady-state  operating  condition  for  fuel  processor  system 


T  (°C) 

CH4  (mol  min  !) 

H20(molmin  !) 

CO  (mol  min  1 ) 

CO2  (mol  min  !) 

O2  (mol  min  !) 

H2  (mol  min  !) 

N2  (mol  min  ]) 

ATR 

Inlet 

111 

0.480 

0.696 

0 

0 

0.216 

0 

0.812 

Outlet 

650 

0.015 

0.413 

0.274 

0.191 

0.03 

1.21 

0.812 

Cooling 

— 

— 

0.486 

— 

— 

— 

— 

— 

HTS1 

Inlet 

350 

0.015 

0.898 

0.274 

0.191 

0.03 

1.210 

0.812 

Outlet 

402 

0.015 

0.710 

0.081 

0.390 

0.03 

1.410 

0.812 

Cooling 

— 

— 

0.200 

— 

— 

— 

— 

— 

HTS2 

Inlet 

317 

0.015 

0.910 

0.081 

0.390 

0.03 

1.410 

0.812 

Outlet 

322 

0.015 

0.880 

0.059 

0.410 

0.03 

1.427 

0.812 

Cooling 

— 

— 

0.220 

— 

— 

— 

— 

— 

LTS 

Inlet 

237 

0.015 

1.110 

0.059 

0.410 

0.03 

1.427 

0.812 

Outlet 

241 

0.015 

1.090 

0.043 

0.420 
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Fig.  3  shows  a  cold  start-up  of  ATR  pathway.  Initially,  the 
reactor  beds  are  at  room  temperature  (25  °C)  and  the  inlet 
of  the  reformer  is  heated  up  with  the  external  fuel  followed 
by  the  recycled  hydrogen.  The  temperature  front  propagates 
through  the  reformer  toward  the  HTS1  after  5  min.  into  the 
start-up.  After  t=  15  min,  the  temperature  front  passes  through 
HTS2  and  to  LTS  at  t-  25  min.  Thirty  minutes  after  the  start¬ 
up,  the  temperature  front  propagates  through  PROX  as  shown 
in  Fig.  3A.  The  total  CO  flow  propagates  at  a  slightly  slower 
rate.  It  reaches  an  almost  static  profile  30  min  after  the  start¬ 
up.  It  is  also  important  to  note  that  the  water  gas  shift  reactors 
(HTS1,  HTS2,  and  LTS)  will  become  effective  when  the  reac¬ 
tor  temperatures  approach  nominal  steady-state  value  (Fig.  3B). 
It  is  interesting  to  note  that  the  CO  concentration  at  the  LTS 
effluent  goes  through  a  non-monotonic  change.  At  15  min, 
the  total  CO  flow  rate  approaches  a  maximum,  because  of 
the  production  of  the  CO  from  the  reformer.  As  the  HST1, 
HTS2,  and  LTS  become  effective,  the  CO  concentration  is 
brought  down  and,  eventually,  reaches  acceptable  level  after 
t>  30  min. 

Note  that  the  temperature  and  concentration  profiles  in  Fig.  3 
is  plotted  in  the  reactor  axial  position  with  the  actual  dimension 
in  the  reactor  length.  Fig.  3  reveals  the  inherent  dynamics  of 
the  fuel  processor  and  it  indicates  without  any  adjustment  in 
operating  condition,  it  takes  more  than  30  min  for  the  reactors 
to  reach  steady- state  profiles. 

3.  Sensitivity  analysis 

Simulation  results  in  Fig.  3  reveal  that  two  factors  play  an 
important  role  in  the  start-up  time  of  a  fuel  processor.  One 
is  the  temperature  of  the  catalyst  bed,  and  the  other  is  the 
composition,  CO  concentration  in  particular,  at  the  effluent  of 
the  PROX  reactor.  Because  the  reactors  are  in  series  cascade 
structure,  if  the  reactor  temperatures  do  not  reach  the  desig¬ 
nated  reaction  temperatures,  the  gas  cleaning  units  may  not  be 
functioning  fully  and  this  results  in  higher  CO  concentration 
at  the  reactor  outlet.  Note  that,  the  fuel  cell  system  will  not 
be  in  operation  if  the  CO  concentration  at  the  outlet  of  the 
PROX  reactor  exceeds  the  specification,  45  ppm.  The  temper¬ 
ature  propagation  in  the  fuel  processor  starts  from  the  ATR. 
Recall  that  we  use  a  furnace  to  preheat  the  fresh  feeds  and  when 
the  temperature  of  the  feeds  reach  the  ATR  reaction  tempera¬ 
ture  (around  700  °C),  hydrogen  is  generated  along  with  CH4, 
H2O,  H2,  CO,  CO2  and  air.  The  hot  outlet  gaseous  stream 
of  the  ATR  serves  as  a  heat  source  which  heats  up  the  cat¬ 
alyst  bed  of  the  following  reactors,  HTS1,  HTS2,  LTS,  and 
PROX. 

Thus,  it  is  important  to  devise  a  strategy  such  that  the  temper¬ 
ature  front  propagates  rapidly  down  the  gas  cleaning  units,  while 
maintaining  low  CO  concentration.  In  other  words,  the  effect  of 
CO  concentration  should  also  be  taken  into  consideration  when 
one  is  devising  a  workable  start-up  strategy.  Without  chang¬ 
ing  the  process  configuration,  available  manipulated  variables 
include:  CH4  feed  flow,  steam  feed  flow,  and  air  feed  rate.  By 
manipulated  variables,  we  mean  the  variables  can  be  adjusted, 
typically  via  control  signal  applied  on  a  control  valve.  They 


can  be  flow  rates  or  ratios  of  two  flow  rates.  The  importance  of 
manipulated  variables  for  rapid  start-up  is  explored  via  sensitiv¬ 
ity  analysis. 

3.1.  Steam-to-carbon  ratio 

First,  let  us  examine  the  manipulated  variable  is  the  “steam- 
to-carbon”  (S/C)  ratio.  The  S/C  is  a  flow  ratio  which  is  adjusted 
by  varying  the  set  point  of  a  ratio  controller  (RC1  in  Fig.  2). 
When  CH4  flow  rate  is  kept  constant,  a  change  in  S/C  implies 
adjusting  the  steam  flow  rate.  Table  1  shows  that  the  steam 
reforming  is  highly  endothermic  and  water  gas  shift  reaction 
is  a  mildly  exothermic  reaction.  An  increase  in  the  S/C  ratio 
favors  the  steam  reforming  reaction  with  less  CO  generation. 
Because  the  inlet  temperature  to  the  reformer  is  controlled  using 
the  furnace  heat  input,  the  heat  (Q  in  Fig.  4)  propagates  down 
the  gas  cleaning  units  will  increase.  That  is  the  heat  capacity 
flow  at  the  outlet  of  the  reformer  increases,  as  a  result  of  energy 
content  carried  by  additional  steam.  The  increased  heat  input  to 
the  furnace  can  be  compensated  using  the  recycled  hydrogen  at 
the  start-up  stage.  Also  note  that  the  outlet  temperature  of  the 
reformer  actually  decreases.  On  the  other  hand,  the  CO  concen¬ 
tration  also  decreases  as  a  result  of  lowered  reactor  temperatures 
(Fig.  4).  This  is  important,  because  an  increase  in  the  S/C  ratio 
can  provide  higher  heat  flux  for  the  ATR  outlet  stream,  while 
keeping  the  CO  concentration  down. 

3.2.  Air-to-carbon  ratio 

The  second  likely  manipulated  variable  is  the  air-to-carbon 
ratio  where  the  air  flow  rate  is  adjusted  while  fixing  CH4  and 
steam  flow  rates.  This  manipulated  input  is  used  by  Springmann 
et  al.  [15]  to  speed  up  the  start-up.  Large  excess  of  oxygen  favors 
the  partial  oxidation  reaction  (Table  1)  and  the  heat  flux  toward 
the  gas  cleaning  unit  goes  up  accordingly.  The  temperature  at  the 
ATR  outlet  also  goes  up  as  a  result  of  highly  exothermic  partial 
oxidation  reaction  as  shown  in  Fig.  4.  However,  a  higher  reactor 
temperature  results  in  a  higher  equilibrium  conversion  for  the 
steam  reforming  (generating  more  CO)  and  a  lower  conversion 
for  the  exothermic  water  gas  shift  reaction  (converting  less  CO 
to  CO2).  Thus,  we  have  a  higher  CO  concentration  as  shown 
in  Fig.  4.  Unlike  changing  S/C  ratio,  adjusting  the  air-to-carbon 
ratio  has  a  positive  effect  on  the  heat  flux  ( Q )  and  a  less  favorable 
effect  on  CO  concentration. 

3.3.  CH4flow  rate 

Let  us  first  examine  the  case  when  the  methane  feed  flow  rate 
is  varied  while  keeping  the  water  and  air  flow  rates  fixed.  Similar 
to  the  case  of  adjusting  air  flow  rate,  both  the  heat  flux  (g)  and 
ATR  outlet  CO  concentration  increase  (Fig.  4).  The  reasons  for 
that  is  a  higher  CH4  concentration  favors  the  steam  reforming 
and  partial  oxidation  reaction  (Table  1),  so  the  temperature  at  the 
ATR  outlet  increases  and  the  CO  generated  also  shows  a  trend 
of  increase. 

An  alternative  is  changing  the  CH4  flow  rate  while  keep¬ 
ing  the  S/C  and  air-to-carbon  ratios  constant.  Certainly,  the  heat 
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Fig.  4.  Sensitivities  of  heat  capacity  flow  rate  ( Q )  and  CO  concentration  at  the  outlet  of  the  ATR  reactor  by  changing  (A)  air-to-carbon  ratio,  (B)  CH4  flow  rate,  (C) 
steam-to-carbon  ratio,  (D)  CH4  flow  rate  while  fixing  air-to-carbon  and  steam-to-carbon  ratios. 


flux  ( Q )  increases  as  a  result  of  an  increase  in  the  throughput. 
But  the  CO  concentration  is  kept  rather  flat  as  the  methane 
flow  rate  changes  (Fig.  4).  Another  advantage  of  using  this 
approach  is  that  the  feed  ratios  of  start-up  are  the  same  as 
that  of  the  normal  operation.  Thus,  when  switching  back  to 
the  mode  of  nominal  steady-state,  only  flow  rate  disturbance 
is  experienced,  instead  of  feed  rate  and  feed  ratio  disturbances. 
Because  we  typically  design  the  nominal  operation  condition 
with  the  flow  rate  set  to  a  certain  percentage  (e.g.,  50%)  of  the 
maximum  flow  rate,  keeping  constant  feed  ratios  is  less  likely 
prone  to  constraint  violation  while  utilizing  the  system  capacity 
fully. 


4.  Start-up  strategies 

4.1.  Control  structure 

Before  getting  into  the  start-up  strategies,  the  control  struc¬ 
ture  for  normal  operation  is  examined  first.  Lin  et  al.  [17] 
proposed  two  control  structures  to  handle  hydrogen  production 
rate  variations.  One  is  the  “on- supply”  control  structure  where 
the  CH4  feed  flow  is  adjusted  directly  to  meet  the  load  demand. 
The  other  is  the  ”on-demand“control  structure  [20,21]  where  the 
hydrogen  flow  rate  at  the  exit  of  the  PROX  reactor  is  set  by  the 
fuel  cell  demand  and  a  pressure  control  is  used  to  set  the  CH4 
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Fig.  5.  Transient  responses  of  inlet  and  outlet  temperatures  and  component  flow  rate  of  gas  cleaning  unit  using:  nominal  feedback  structure  (dashed)  and  optimized 
feed  manipulation  strategy  (solid). 


feed  flow  rate.  The  same  control  structure  is  used  here,  except 
that  the  pressure  controller  is  put  in  manual  mode  at  the  start-up 
stage.  Taking  out  a  flow  loop  and  a  pressure  loop,  the  control 
structure,  at  the  start-up  stage,  consists  of  the  following  loops 
(Fig.  2): 

1.  CH4  feed  is  under  flow  control  (FC  in  CH4  flow  path  in 
Fig.  2). 

2.  H2O  flow  rate  is  ratio  to  CH4  feed  rate  (RC1  in  Fig.  2). 

3.  O2  flow  rate  is  ratio  to  the  CH4  feed  flow  (RC2  in  Fig.  2). 

4.  Reformer  inlet  temperature  is  maintained  by  changing  the 
fuel  flow  rate  (TC  providing  the  set  point  of  the  flow  con¬ 
troller). 

5 .  HTS 1  (first  high-temperature  water-gas  shift  reactor),  HTS2, 
LTS  inlet  temperatures  are  controlled  by  water  injection. 

6.  PROX  inlet  temperature  is  maintained  by  varying  water  injec¬ 
tion  rate. 

Fig.  2  shows  the  control  structure  with  three  flow  loops,  two 
ratio  controllers,  one  temperature  loop  for  fuel  flow,  four  tem¬ 
perature  loops  for  H2O  addition.  Note  that  instrument  signals 
are  depicted  is  dashed  line  in  Fig.  2.  At  the  start-up  stage,  the 
composition  loop  does  not  work  until  the  process  reaches  to 
steady-state.  Here,  the  set  point  of  RC1  and  RC2  are  fixed.  A 
third  order  lag  with  a  time  constant  of  0.5  min  is  assumed  for  the 
reformer  inlet  temperature  and  the  relay  feedback  test  is  used  to 
find  the  ultimate  gain  and  ultimate  period  [22].  The  PI  controller 
settings  are  obtained  from  Tyreus-Luyben  tuning  method  [23]. 
In  order  to  provide  a  basis  for  comparison,  let  us  first  exam¬ 
ine  the  start-up  time  using  the  nominal  control  structure.  That 


is  all  controllers  are  set  to  automatic  when  the  fuel  processor  is 
switched  on,  i.e.,  without  any  designated  start-up  strategy.  Simu¬ 
lation  results,  dashed  line  in  Fig.  5,  shows  that  it  takes  37.5  min 
to  reach  operable  steady-state,  i.e.,  CO  concentration  meeting 
the  45  ppm  spec.  Note  that  the  long  start-up  time  is  due  to  large 
thermal  masses  of  the  HTS  1 ,  HTS2,  LTS  and  PROX  reactors  and 
the  limitation  on  the  capacity  of  fuel  flow  of  the  experimental 
system.  Despite  all  these  limitations,  it  is  possible  to  speed-up 
the  cold  start-up  by  feed  manipulation. 

4.2.  Feed  manipulation 

Sensitivity  analysis  in  Section  3  reveals  two  important  factors 
to  reduce  the  start-up  time:  (1)  making  the  heat  flux  Q  as  large  as 
possible  and  (2)  providing  smooth  transition  back  to  the  normal 
operating  point.  Looking  at  the  steady-state  heat  flux  and  CO 
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Fig.  6.  Optimization  variables,  pulse  height  ( h ),  pulse  duration  (tf),  and  ramp 
down  slope  ( S ),  for  start-up  time  minimization. 
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concentration  effects  in  Fig.  4,  the  CH4  feed  flow  rate  with  fixed 
S/C  and  air-to-carbon  ratios  seems  to  be  an  ideal  candidate  for 
reducing  the  start-up  time.  Assuming  a  pulse  change  in  the  feed 
flow  as  shown  in  Fig.  6,  then  the  questions  becomes: 

1.  What  is  the  pulse  height  (h)l 

2.  What  should  be  the  pulse  duration  time  (/a)? 

3 .  What  is  the  appropriate  slope  ( S )  ramping  down  to  the  normal 

operating  point? 

The  answer  to  the  first  question  is  quite  obvious.  In  order  to 
get  as  much  heat  flux  as  possible,  the  feed  is  set  to  its  maximum 
value,  e.g.,  FCh4  =  TcH4,max.  In  a  typical  engineering  design, 
the  control  valve  is  half-open  at  nominal  flow  rate.  So  the  feed 
flow  is  doubled  (from  its  nominal  value)  initially  (Fig.  6). 

Once  the  pulse  height  ( h )  is  fixed,  the  effect  of  the  duration 
time  (/a)  is  explored  next  (Fig.  6).  Let  us  use  square  pulses, 
i.e.,  ramp  down  slope  S  approaching  infinity,  to  illustrate  the 
effects  of  duration  time  to  the  cold  start-up.  Fig.  7  shows  dynamic 
responses  temperatures,  flows  and  concentration  for  three  dif¬ 
ferent  duration  times.  As  the  pulse  duration  time  increases,  the 
temperature  front  propagates  down  the  reactors  (e.g.,  outlet  tem¬ 
perature  of  PROX  reactor)  at  a  faster  speed.  But  a  negative 
step  change,  back  to  the  nominal  operating  point,  introduces 
additional  disturbance  which  subsequently  driving  the  CO  con¬ 


centration  out  of  spec.  Fig.  7  shows  that  for  the  duration  time 
of  10,  18,  and  35  min,  the  corresponding  start-up  times  are 
41.3,  29.1,  and  33.7  min,  respectively.  Two  observations  can  be 
made  immediately.  First,  a  negative  step  change  (resulted  from  a 
square  pulse)  is  too  large  a  disturbance  for  a  stabilizing  system. 
Second,  there  exists  an  optimal  duration  time.  That  is  too  small  a 
duration  time  will  slow  down  the  temperature  front  propagation 
and  too  large  a  duration  time  makes  the  start-up  time  effective 
the  same  as  the  start-up  with  the  nominal  control  structure  (e.g., 
dashed  line  in  Fig.  5). 

Thus,  the  effects  of  ramp  down  slope  are  investigated. 
Assume  the  duration  time  is  t^  =  18  min,  three  different  ramp 
down  rates,  (S=  AF/At=  —  00  (—0.48/0),  —0.096  (—0.48/5), 
and  —0.024  (—0.48/20)),  are  simulated  as  shown  in  Fig.  8.  Fig.  8 
shows  that  for  the  duration  slope  of  00,  —0.096,  and  —0.024  min, 
the  corresponding  start-up  times  are  29.1,  28.8,  and  29.8  min, 
respectively.  The  results  indicate  that  a  large  a  slope  introduces 
unnecessary  disturbance  in  CO  concentration  and  too  small  a 
slope  lead  to  a  sluggish  response.  Again,  an  optimal  setting  exists 
for  the  ramp  down  rate. 

Instead  of  going  through  all  possible  combinations,  the  start¬ 
up  problem  can  be  formulated  as  an  optimization  problem, 
provided  the  shape  of  feed  manipulation  (e.g.,  Fig.  6).  The  objec¬ 
tive  is  to  minimize  the  start-up  time  (/startup)  which  is  defined  as 
the  time  where  CO  concentration  at  PROX  outlet  reaches  45  ppm 
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Fig.  7.  Cold  start-up  of  the  fuel  processor  with  pulse  duration  time  of  10  min  (dashed),  18  min  (solid),  and  35  min  (dot). 
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Fig.  8.  Cold  start-up  of  the  fuel  processor  with  ramp  down  slope  of  -oo  (dot),  —0.096  (solid),  and  —0.024  (dashed). 


and  remains  below  45  ppm  for  t  >  /startup-  The  optimization  prob¬ 
lem  can  be  expressed  as 

iftin /startup  —  f(h,  /fu  3),  £2  =  {h,  S} 

subject  to  {h  <  /zmax} 

where  h  is  the  pulse  height,  t&  the  duration  time,  S  the  ramp  down 
slope,  and  /zmax  is  the  allowable  input  change.  The  results  show 
that  optimal  t&  and  S  are  17  min  (to)  and  —0.096  (S'),  respec¬ 
tively.  With  this  setting,  simulation  result,  solid  line  in  Fig.  5, 
shows  that  the  start-up  time  can  be  reduced  by  a  factor  of  25% 
(as  compared  to  the  nominal  control  structure-based  start-up, 
/startup  =  37.5  min)  and  becomes  /startup  =  28  min. 

4.3.  Heuristics 

Despite  the  fact  that  mathematical  programming  can  be  used 
to  find  the  optimal  start-up  strategy,  it  will  be  helpful  if  an 
efficient  start-up  can  be  generated  heuristically.  Extensive  sim¬ 
ulation  reveals  that  a  suboptimal  solution  can  be  obtained  using 
the  following  approach.  Provided  with  the  nominal  operating 
condition,  we  have  the  following  steps: 

1.  Set  the  CH4  feed  flow  to  the  maximum  value.  Note  that  the 
air  and  water  flow  rates  also  increase  with  the  same  ratio. 


2.  Identify  the  time  (i.e.,  time  =  framp)  when  the  PROX  inlet 
temperature  starts  to  increase  and  also  find  the  rate  of  the 
temperature  change  (i.e.,  A7prox/A0. 

3 .  Start  to  ramp  down  at  framp  with  the  slope  A  Fcu4  /  [(?prox  — 
25)/2/(A7prox/A0]-  That  is  the  ramp  down  is  completed 
when  the  reactor  inlet  temperature  is  half-way  (i.e.,  (7prox  — 
25)/2)  to  the  steady-state  value  (7prox)- 

With  this  start-up  strategy,  the  result  shows  that  the  start-up 
dynamics  is  almost  the  same  as  that  of  the  optimal  one  (solid 
line  in  Fig.  5)  and  the  start-up  time  is  28  min. 

4.4.  Design  modifications 

Up  to  this  point,  we  have  devised  start-up  strategies  for  an 
existing  fuel  processor  (Fig.  1).  The  next  question  then  becomes: 
can  we  speed  up  the  start-up  time  further  by  design  modification. 
From  design  perspective,  two  approaches  are  possible.  One  is 
to  use  the  heat  from  the  effluent  of  PROX  reactor  to  heat  up  the 
inlet  streams  to  HTS1,  HTS2,  LTS,  and  PROX,  a  typical  heat- 
integration  approach  (Fig.  9).  This  can  speed-up  the  temperature 
front  propagation  by  providing  heat  at  the  reactor  inlet.  The  sec¬ 
ond  approach  is  associated  with  materials  design.  Recall  that  the 
inherent  problem  with  large  start-up  time  of  the  fuel  processor  is 
the  heating  rate.  Because  of  the  heat  capacity  of  the  air  is  quite 


PROX.irr  *  LTS,inv  '  HTS2,irr  '  HTS1,in 
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small,  the  convective  heating  dominates  the  start-up  time.  If  low 
heat  capacity  (\ow-Cp)  materials  can  be  used  for  the  catalyst  car¬ 
rier,  the  start-up  time  can  be  reduced  significantly.  Consider  the 
case  when  the  heat  capacity  of  the  carrier  is  reduced  to  50%  of 
their  original  values.  Simulation  results,  Fig.  10,  show  that  tem¬ 
perature  front  is  propagated  at  a  much  faster  rate  and  the  start-up 
time  becomes  22  min,  a  42%  reduction  from  the  original  design. 
Note  that  this  is  achieved  using  the  nominal  control  structure 
for  the  start-up  and  further  improvement  is  possible  using  the 
optimal  feed  policy  as  discussed  earlier. 

5.  Conclusion 

Cold  start-up  of  a  methane  fuel  processor  is  studied  in  this 
work.  The  experimental  fuel  processor  is  intended  to  provide 
hydrogen  for  a  proton  exchange  membrane  (PEM)  fuel  cell  for 
the  power  generation  (3  kWe).  A  dynamic  model  consisting  of 
a  reformer,  three  water-gas  shift  reactors,  and  a  preferential 
reactor  are  described.  Two  determining  factors  for  a  rapid  start¬ 
up  are  identified:  speed  of  temperature  front  propagation  and 
acceptable  CO  concentration.  Steady-state  analyses  reveal  that 
the  fuel  feed  flow  rate  with  fixed  steam-to-carbon  and  air-to- 
carbon  ratios  is  an  ideal  manipulated  variable  at  the  start-up 
stage.  Considering  both  large  initial  heat  flux  for  convective 
heat  transfer  and  a  gradual  transition  back  to  nominal  operat¬ 
ing  point,  the  shape  of  feed  manipulation  is  determined  which 
is  a  pulse  with  a  linear  ramp  down  to  the  nominal  value.  With 
the  feed  scenario  available,  the  fuel  processor  start-up  can  be 
formulated  as  a  constrained  optimization  problem  and  can  be 
solved  numerically.  From  optimization  result,  a  heuristic  is  gen¬ 
erated  for  rapid  start-up  of  a  fuel  processor.  This  leads  to  a  25% 
improvement  in  the  start-up  time.  Finally,  issues  of  design  mod¬ 
ification  are  explored  for  further  reduction  in  the  start-up  time. 
The  results  show  that  if  the  low  heat  capacity  materials  are  used 
for  catalyst  carriers,  the  start-up  time  can  be  reduced  by  a  factor 
of  42%. 
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